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1.  INTRODUCTION 

This  second  semi-annual  report  summarizes  die  fHogress  made  in  the  last  six  months.  The  overall 
objectives  of  diis  research  are  to  develop  new  methods  to  evaluate  subgrid  models  and  then  to  uti¬ 
lize  these  methods  to  improve  the  chosen  subgiid  models.  The  subgiid  models  investigated  in  this 
research  are  chosen  primarily  for  triplication  in  high  Reynolds  number  complex  flows.  Prelimi¬ 
nary  studies  of  these  models  have  b^  completed.  A  priori  analysis  u;>mg  data  from  direct  numer¬ 
ical  simulations  (DNS)  homogeneous  isotropic  flows  was  carried  out,  and  tl^n  the  models  were 
imidemented  in  large-eddy  simulation  (LES)  codes  and  fuidier  evaluated.  IWo  types  of  analysis 
mefliods  have  been  deveioped  so  far.  The  first  method  uses  information  in  Fourier  (qiectral)  qiace 
and  evaluates  the  interscale  energy  tranter  as  a  function  of  the  wavenumbers  resolved  in  the  LES. 
Ibe  sectmd,  uses  infoimation  in  flie  {diysical  space  and  uses  cross  correlation  analysis  to  investi- 
gate  the  behavior  of  subgrid  models.  Tte  physical  ^ace  analysis  method  will  be  the  primary  anal¬ 
ysis  tool  for  the  next  year’s  study,  since  ^  next  ph^  of  research  will  focus  on  complex  flows 
such  as  flows  past  rearward  facing  steps  and  swiriing  flows. 

bi  the  following,  we  discuss  recent  devdopments  of  the  last  six  months.  Some  of  the  results  have 
been  described  in  more  detail  in  the  ptqrers  (Moron  and  Yeung,  1994a,  1994b;  Kim  et  al.,  1994) 
and,  therefore,  will  not  be  rqreated  here.  These  papos  (or  the  extended  abstracts)  are  included  as 
ai^rordices  to  this  report 

2.  ANALYSIS  OF  SUBGRID  MODELS  USING  DNS  AND  LES  OF  ISOTROPIC  TURBU- 
UBMCE 

bi  dris  phase  of  research,  homogoieous  isotropic  turbulence  was  used  as  a  test  flow  field  to 
devdop  analysis  methods.  This  was  motivated  by  the  fact  that  isotropic  turbulence  has  been  stud¬ 
ied  extoisivdy  in  the  past  and  detailed  DNS  data  is  available.  In  addition,  the  computational 
domain  is  very  simple  and  ttoefore,  medrods  using  qrectral  techniques  can  be  used  to  identify 
features  of  the  flow  which  are  not  possible  using  {diysical  qrace  methods  and  vice  versa.  Both 
hicompresrible  and  compressible  isotropic  fields  have  been  studied  using  two  different  finite  dif¬ 
ference  codes.  These  cotfes  have  been  validated  by  comparing  die  DNS  predictions  by  these  codes 
widi  the  results  obtained  using  the  well  know  qiectral  code  of  Rogallo. 


2.1  Analysis  of  Incompressible  Isotropic  Turbulence 

^Auioos  subgrid  models  were  evaluated  in  both  spectral  and  physical  ^ace  using  high  resolution 
DNS  data.  Subsequently,  these  models  were  implemented  in  coarse  grid  LES.  Tbe  ^)ectral  ^ace 
mediod  was  used  primi^y  for  a  priori  analysis  while  the  physical  space  method  was  used  for 
badsapriori  and  posteriori  analysis.  For  a  priori  analysis,  vesy  high  resolution  data  (using  128^ 
grid)  was  used.  All  the  posteriori  analysis  were  carried  out  using  LES  data  and  without  using  any 
DNS  infoimatitm.  This  iq^Hroach  isessmitial  fen*  future  application  to  high  Reynolds  number  flows 
since  DNS  data  for  such  flows  will  not  be  available. 

For  the  physical  q;>ace  analysis  of  the  subgrid  models,  LES  using  different  grid  resolutions  is  first 
carried  out  Then,  using  top  hat  filtering,  the  modeled  subgrid  stresses  and  the  energy  transfer  are 
correlated  betwem  the  two  LES  data  fields.  For  example,  LES  was  carried  out  using  32^  grid 
and  16^  grid  using  identical  initial  flow  field  and  using  the  same  subgtid  model.  Then,  at  a  chosen 
instant,  the  resolved  field  in  the  32^  grid  is  filtered  to  compute  the  subgrid  stresses  and  energy 
transfer  in  die  16^  grid.  This  field  is  considered  ‘exact',  as  far  as  the  coarser  grid  is  concerned, 
since  the  length  scales  that  are  unresolved  (and  hmice  modeled)  in  fire  coarser  grid  are  supposed  to 
be  resolved  in  the  finer  grid  LES.  Therefore,  if  the  model  is  wmlring  properly  this  ‘exact’  field 
must  be  tqrroduced  by  the  subgrid  model  in  the  16x16x16  grid.  Cross  correlation  analysis 
between  die  modeled  and  ‘exact  subgrid  stresses  and  the  energy  transfer  was  carried  out  If  the 
cmrelation  is  high,  it  would  suggest  fiiat  the  subgrid  model  behaves  ctmsistently  for  different  grid 
resolutions  and  that  die  subgrid  energy  transfer  is  modeled  correedy  in  the  16x16x16  LES.  With 
this  iqiproach,  model  validation  does  not  require  DNS  data  and,  more  importandy,  since  the  same 
flow  fidd  is  being  investigated,  die  model  can  be  investigated  diiecdy  in  the  flow  field  and  geom¬ 
etry  of  interest  Furthermore,  this  qiproach  allows  an  immediate  assessment  of  the  capability  of 
die  subgrid  models  in  the  coarsest  grid. 

Note  that  die  results  of  the  above  analysis  methods  do  not  provide  any  information  on  the  accu¬ 
racy  of  the  results.  Comparison  with  experimental  data  (or  DNS  data,  uiiere  ever  possible)  is 
essential  to  demonstrate  the  accuracy  of  the  LES.  So  far,  for  isotix^c  flows,  DNS  data  have  been 
used  to  evaluate  the  accuracy  of  the  L£S  results,  but  future  studies  will  be  directed  to  more  com- 
idex  flows  for  which  no  DNS  data  is  availaUe.  Cmnparison  widi  experimental  data  will  have  to 
be  carried  mit  in  such  cases.  It  is  expected  that  the  analysis  methods  fm  model  validation  will  also 
have  to  be  further  developed  for  cennplex  flows,  fenr  example,  to  handle  the  near  wall  effects. 

Vhrious  subgrid  models  have  been  implemmited  and  evaluated  using  the  techniques  described 
above.  More  details  of  the  analysis  methods  ate  given  in  Mmion  and  Yeung  (1994a,  1994b).  The 
subgrid  models  studied  so  fxr  are: 

(a)  die  classical  Smagminsky's  eddy  viscosity  model 
^)  die  dynamic  (Germano's)  eddy  viscosity  model 

(c)  a  qiectral  edtfy  viscosity  model 

(d)  a  new  scale  similarity  model 

(e)  a  one-equadtm  model  for  the  subgrid  kinetic  energy  with  and  without  stochastic  badcscatter 

(f)  a  dynamic  tme-equation  model  for  the  subgrid  kinetic  energy 

^)  a  two-equation  model  for  die  subgrid  kiiietic  energy  and  subgrid  helicity  (k-h  model) 


Smagorinsky’s  model  is  voy  popular  in  literature;  however,  it  has  been  ^own  to  require  signifi¬ 
cant  modificatioas  (primarily  adjustment  of  the  ‘constant’)  for  good  agreement  with  experimental 
data.  Hie  major  ‘breaktinough’  in  subgrid  model  development  is  the  application  of  the  algebraic 
identity  of  Gomano  to  evaluate  dynamically  die  ctmstant  in  Smagorin^’s  model.  In  Menon  and 
Yeung  (1994a),  we  studied  the  classical  mo^l  widt  fixed  constant  while  in  Menon  and  Yeung 
(1994b)  and  Kim  et  al.  (1994),  we  studied  the  dynamic  eddy  viscosity  model. 

Ihe  scale  similarity  model  is  a  modified  vcasion  of  the  original  Bardina’s  model  and  was  pro¬ 
posed  by  Meneveau  (Liu  et  al.,  1993)  based  on  analysis  of  high  Reynolds  number  experimental 
data  on  turbulent  jets.  It  is  based  on  the  idea  that  the  energy  transfer  at  the  resolved  grid  resolution 
is  self  similar  to  the  energy  transfer  occuring  at  a  resolution  twice  as  coarse.  This  method,  there- 
fme,  uses  a  test  filter  (an  apimach  very  similar  to  Germano’s)  to  compute  the  scale-similar  sub¬ 
grid  stresses  in  terms  of  the  resolved  field.  Computationally,  t^  model  is  very  simple  and  easy  to 
implement  However,  as  discussed  in  Menon  and  Yeung  (1994a).  there  are  some  inhenant  limita¬ 
tions  to  diis  model  when  used  in  LES.  This  model  can  predict  backscatter  but  the  amount  of  back- 
scatter  may  exceed  the  real  backscatter.  This  can  result  in  numerical  instability  and,  therefore, 
smne  smt  of  backscatter  control  is  necessary.  More  details  of  the  analysis  of  this  model  are  given 
below  and  in  Menon  and  Yeung.  (1994a). 

The  (»e  equaticHi  model  for  the  subgrid  kinetic  miergy  (k-equation  model)  was  chosen  keeping  in 
mind  the  requirements  for  practical  high  Reynolds  number  LES.  It  is  expected  that  for  high  Rey¬ 
nolds  number  LES  of  complex  flows,  the  grid  resolution  practically  possible  (due  to  resource  con¬ 
straints)  will  be  limited.  Therefore,  simple  dissipative  models  (even  with  dynamic  evaluations) 
may  not  be  sufficient  for  practical  LES.  In  addition,  die  assumpdcm  of  local  equilibrium  between 
die  production  and  dissipation  of  the  kinetic  energy  (an  assumption  implicit  in  all  algebraic  eddy 
viscosity  models)  may  be  violated.  The  k-equation  model  with  fixed  coefficients  was  investigated 
in  Mentm  and  Yeung  (1994a),  while  the  dynamic  k-equation  model  is  investigated  in  Menon  and 
Yeung  (1994b)  and  Kim  et  al.  (1994). 

We  are  also  investigating  more  advanced  models  for  high  Reynolds  number  flows.  One  such 
model  is  a  two-equation  model  for  the  subgrid  kiiretic  energy  and  subgrid  helicity  (the  k-h  mottel). 
Helicity  is  non-zero  only  if  die  flow  is  locally  3D.  Thus,  if  the  small  scales  are  anisotropic  or  non- 
homogeneous  (which  can  occur  if  die  grid  is  coarse,  the  geometry  is  complex,  and  the  Reynolds 
numbm’  is  very  high)  then  simple  eddy  viscosity  models  or  even  the  one-equation  models,  may 
not  be  able  to  take  into  account  diis  small-scale,  local  3D  effects.  Some  preliminary  studies  of  this 
model  have  been  completed  and  results  are  discussed  below. 

•  ZJJ  Summary  of  the  Results 

The  results  of  die  arudysis  of  diese  sobgrid  models  have  been  rqiorted  in  the  p^rs  attached  in 
the  Appendices.  Here,  we  briefly  summarize  those  results  and  then  discuss  some  new  results 
tecendy  obtained. 

The  analysis  of  the  eddy  viscosity  models  (models  (a)  and  (c)),  the  scale  similarity  model  (model 
(d))  and  die  one-equatkm  model  (model  (e))  ytere  rqiorted  in  Menon  and  Yeung  (1994a).  The  a 
priori  analysis  showed  that  for  fine  grid  resolution,  the  scale  similarity  model  had  the  highest  cor- 


idaticai  with  the  exact  subgrid  stresses  and  energy  transfer.  However,  for  a  coarse  grid  LES,  this 
CMTelatitm  drof^ped  significantly  indicating  that  this  model  is  not  iqppropriate  for  coarse  grid  LES. 
This  result  was  understandable  since  the  scale  similarity  concept  implies  that  the  energy  transfer 
at  two  grid  levels  are  self  similar.  As  die  grid  is  coarsened,  this  similarity  begins  to  breakdown. 
For  die  low  Reynolds  number  flows  studied  in  Mentm  and  Yeung  (1994a).  there  was  no  clear  iner¬ 
tial  range  resolved  in  die  DNS.  This  made  it  difScolt  to  fully  evaluate  this  model.  This  model  was 
proposed  for  high  Reynolds  number  flows  (based  cm  experimental  data  at  Res310)  udiere  a  dis¬ 
tinct  ineiital  range  existed.  Scale-similarity  assumptions  hold  very  well  in  the  inertial  range  and, 
thus,  this  model  would  be  q;^licable.  However,  this  implies  that  to  use  this  model  an  initial 
range  must  be  resolved  in  the  LES.  This  may  be  an  unacceptable  requirement  in  high  Reynolds 
number  flows  (since,  it  implies  a  very  high  grid  resolution)  and.  thus,  at  (nesmit.  diis  model 
(aldiough  very  elegant)  appears  to  have  limited  use  for  LES  of  complex,  high  Reynolds  number 
flows. 

&nagorinsky’s  model  (with  fixed  constant)  was  quite  poor  when  compared  to  the  other  models. 
The  k-equation  model  (with  fixed  coefficients)  was  better  than  the  eddy  viscosity  mottel  but  had  a 
correlation  lowm'  than  the  scale  similarity  moitel.  Interestingly,  when  the  grid  was  coarsei^d.  the 
scale  rimilarity  model  became  poorly  ctxielated;  however,  the  k-equation  model  did  not  show 
such  a  behavior.  This  suggests  that  the  k-equation  model  has  the  potential  for  modeling  the  sub¬ 
grid  stresses  and  enmgy  transfer  in  coarse  grids. 

Subsequent  to  this  a  priori  analysis,  diese  models  were  implemented  in  a  LES  code  and  simula- 
titms  were  carried  out  using  different  grid  resolutions.  Ib  compare  the  results,  all  simulations 
were  begun  with  nearly  identical  initial  conditions.  Thus,  at  tsO,  aU  the  flow  fields  were  highly 
corielaied.  The  analysis  was  carried  out  at  an  instant  when  die  flow  had  evolved  to  realistic  turbu¬ 
lence.  The  cmtelation  analysis  using  just  the  LES  data  showed  a  completely  different  picture.  AU 
models  showed  very  poor  correlation  (compared  io  a  priori  tests).  Both  the  eddy  viscosity  model 
and  dm  scale  similarity  model  appeared  to  model  die  subgrid  stresses  quite  poorly  compared  to 
die  one-equation  model. 

These  results  clearly  highlight  the  fact  that  subgrid  models  that  iqipear  to  be  quite  good  in  a  priori 
analysis  may  not  be  as  good  when  implemtmted  in  actual  LES.  Furthermore,  it  appeared  that  the 
k-equation  model  had  the  best  c^iability  to  model  the  subgrid  stresses  in  coarse  grid  LES. 

Ihe  stiMly  in  Menon  and  Yeung  (1994a),  employed  the  models  with  fixed  coefficients.  Since  it  had 
beenshownby  odier  researchers  that  the  dyna^  eddy  viscosity  model  is  quite  superior  to  Sma- 
gofinsky’s  clasrical  model,  it  was  decided  to  revisit  tb^  models  but  by  aUowing  for  dynamic 
evalutaiion  of  the  co^cients.  The  dynamic  eddy  viscosity  model  and  a  dyruunic  version  of  the  k- 
equation  model  was  dien  analyzed  using  the  methods  develt^ied  above.  Some  of  these  results  wiU 
be  reported  in  details  in  Menon  and  Yeung  (1994b)  and  Kim  et  al  (1994).  It  has  been  shown  diat 
die  dyruunic  procedure  significanly  improves  die  correlations.  An  interesting  observation  was 
that,  compart  to  die  dyruunic  Oemano*s  model,  the  dynamic  k-equation  model  showed  a  much 
beaex  improvement  and  was  clearly  supaior  for  coarse  grid  LES.  Ibis  has  givmi  confidence  diat 
for  coarse  grid  LES  fmrhigh  Reynolds  number  flows,  the  use  of  such  higher  order  models  may  be 
beneficial.  This  is  an  issue  that  will  be  revisited  using  more  complex  flows  in  the  second  year  of 
diis  research. 


We  are  now  investigating  a  new  two-equation  model.  Tbe  governing  equations  are  shown  in  the 
Appmidix  and  the  subgrid  model  is  essentiaUy  the  model  prcq)osed  by  Yoshizawa  (1993).  We 
expect  that  in  the  course  of  our  study,  this  model  (if  useful)  will  undmgo  some  changes.  For 
exmni^,  so  far,  results  have  bemi  otKained  using  fixed  coefficients.  We  expect  the  dynamic  pro- 
cedme  will  improve  fliis  model  and,  therefore,  we  are  now  in  the  process  of  including  a  dynamic 
procedure  to  solve  this  model. 

Wb  simulated  a  simple  periodic  flow  field  for  die  ThylOT-Green  Vortex.  This  field  is  (Himarily  2D 
and  will  not  ctmtain  any  helicity.  This  was  confirmed  by  carrying  out  LES  using  the  k-h  model 
and  showing  that  die  subgrid  helicity  was  negligible,  hfext,  the  spanwise  velocity  field  was 
dumged  by  adding  a  cosx  term.  This  allows  the  large-scales  to  becmne  3D  while  sdll  satisfying 
contmuity.  We  were  interested  in  determining  if  this  flow  field  would  generate  small-scale  3D 
structures  and,  if  so,  would  the  new  k-h  model  be  able  to  {mdict  the  helicity  in  the  unresolved 
scales. 

We  lodted  at  three  quantities:  (1)  the  correlation  between  subgrid  helicity  and  large-scale  vortic- 
ity,  (2)  the  accuracy  with  which  subgrid  helicity  is  being  modeled,  and  (3)  the  relevance  of  helic¬ 
ity  to  the  subgrid  stresses  (that  is,  the  effect  of  non-eddy  viscosity  terms  in  the  stress  model  that 
appears  explicidy  due  to  helicity). 

In  Hg.  la  and  lb,  we  show  3D  visualization  of  the  vortex  tubes  (constant  voiticity  isosurfaces) 
along  widi  contours  of  subgrid  helicity  as  predicted  by  the  model.  The  subgrid  helicity  is  only 
shown  on  a  cotain  plane.  Hgure  la  shows  the  field  as  seen  in  a  32x32x32  LES,  while  in  Fig.  lb, 
the  results  are  shown  for  a  16x16x16  LES.  It  can  be  seen  diat  regions  of  intense  vorticity  are  sur¬ 
rounded  by  regions  of  high-magnitude  subgrid  helicity.  This  indicates  significant  production  of 
small-scale  helicity  by  the  breakdown  of  the  large  scale  structures.  This  result  also  implies  that  the 
unresolved  scales  may  be  locally  anisotropic. 

C(»relation  between  the  large-scale  helicity  (due  to  the  resolved  fields)  and  the  subgrid  helicity 
was  also  computed.  This  correlation  should  be  very  low,  in  fact,  it  should  ideally  be  zero,  since 
die  modd  is  siqiposed  to  compute  the  subgrid  helicity  only  due  to  the  anisotropy  or  non-homoge- 
nei^  in  the  small-scales  and,  therefore,  should  not  correlate  with  the  large-scale  helicity.  The 
computed  ccnrelations  were  also  very  small  with  the  32x32x32  LES  predicting  a  value  of  1.03E-4 
and  die  16x16x16  LES  predicting  a  value  of  S.6E'4.  These  results  showed  that  the  k-h  model  has 
been  implmnented  correcdy  and  iqipears  to  be  predicting  die  cmrect  physics. 

Dmrdition  between  die  subgrid  bdid^  modeled  in  the  16x16x16  LES  and  the  subgrid  helicity 
predicted  by  filtering  die  32x32x32  LES  data  into  a  16x16x16  grid  was  also  carried  out  If  the 
model  behaves  accurately  in  both  grids  dien  this  correlation  should  be  high.  Our  fnelimiruuy  data 
diowed  acmielation  of  0.736  which  is  reasonably  high.  Figure  2a  shows  die  subgrid  helicity 
Goaqniied  using  the  32x32x32  LES  data  fihmed  to  16x16x16  grid,  and  Fig.  2b  shows  die  morfel 
prediction  in  the  16x16x16  LES.  dearly,  die  model  is  behaving  reasonable  wdl  in  both  grids. 


This  study  is  not  yet  cmnpleted  and  diere  are  sdll  some  unresolved  issues.  For  example,  the  inclu- 
skm  oi  helicity  model  did  not  improve  the  subgrid  energy  transfer  correlation.  However,  our  pre¬ 
vious  siwty  u^g  fixed  coefficients  for  die  k-equation  model  also  showed  poor  correlation(Menon 


ttul  Yeung,  1994a),  with  dynamic  evaluation,  the  conelation  improved  significantly  (Kim  et 

al.,  1994).  Therefore,  we  are  now  beginning  to  evaluate  dynamically  the  constants  in  this  k-h 
model 


Note  that,  adding  one  more  equation  will  increase  die  cmnputational  cost  Thmfore,  befme  such 
models  are  pn^iosed  for  LES  implication,  it  must  be  cleariy  dmnonstrated  that  it  is  mpeiior  to  die 
conventional  e^  viscosity  mo^.  The  tests  using  Thylor-Green  vmtex  or  isotropic  turbulence 
may  not  be  aipropriate  to  evaluate  this  model  Therefore,  we  are  now  starting  to  implement  this 
model  into  die  cote  developed  to  simulate  more  cmnplex  flows  such  as  flows  past  rearward  facing 
steps.  If  diis  model  is  superior  for  such  flows,  then  die  additional  cost  of  c(»nputation  may  be  bal- 
an^  by  the  ability  of  tte  new  motel  to  handle  amplex  flows  using  relatively  coarser  grids 
(dier^,  decreasing  computational  cost).  This  is  the  primary  goal  of  this  research. 

2.2  Anafysis  of  ComprtssibU  Isotropic  TUrbulence 

Some  studies  were  also  carried  out  to  extend  the  analysis  methods  to  study  compressible  flows. 
As  noted  before,  the  analysis  methods  are  supposed  to  be  intependant  of  the  type  of  flow  field 
studied  and,  therefore,  with  some  minor  mo^cations  should  be  applicable  in  compressible 
flows.  Tb  compare  with  the  incompressible  flow  results,  we  began  by  simulating  low  Mach  num¬ 
ber  (essentially  incompressible)  isotropic  turbulmit  flow  fidds.  So  far,  tmly  the  compressible  ver- 
sicms  of  the  eddy  viscosity  motel  (Eriebacher  et  al,  1987),  the  dynamic  eddy  viscosity  motel, 
and  die  scale  similarity  motel  have  been  implmnented  and  evaluated.  Here,  we  will  summarize 
some  of  the  mme  recent  results  of  this  study.  Mme  details  of  this  work  will  be  included  in  the 
final  version  of  the  paper  Menon  and  Yeung  (1994b). 

Figures  3a  and  3b  show,  respectively,  the  ctnielation  of  tte  exact  subgrid  stress  (obtained 
from  64^  DNS  data)  with  die  eddy  viscosity  motel  and  scale-similarity  model  predictions  as  a 
function  of  filter  width.  Again,  as  before,  box  filters  have  been  employed.  The  results  for  the  ear¬ 
lier  incompressible  data  are  also  shown.  These  figures  show  the  characteristic  decrease  in  correla¬ 
tion  when  die  grid  is  coarsened  with  the  scale  similarity  modd  showing  most  rapid  decrease.  The 
compressible  data  is  quite  dose  to  the  incmnpressible  data  since  very  low  Mach  number  flow  has 
been  simlated.  However,  note  that,  two  completely  difformit  numerical  solvers  and  subgrid  mod¬ 
els  were  employed  for  this  comparison. 

Figures  4a  and  4b  show,  respectivdy,  die  correlation  between  exact  energy  transfer  x.^^J  and  die 
modeled  energy  oansfer  fm'  the  two  models.  Again,  both  models  show  that  with  decrmise  in  grid 
lesolutxm,  die  correlaticm  decreases.  The  scale  similarity  motel  again  shows  a  strong  dependence 
on  the  grid  resoluticm.  However,  for  relatively  fine  mesh,  the  scale  similarity  motel  is  quite  supe- 
lior  to  the  eddy  viscosity  motel.  This  is  in  agreement  with  the  incompressible  flow  results  dis¬ 
cussed  in  Menon  and  Yeung  (1994a). 

ISnce  die  incmnpressible  study  showed  that  the  dynamic  subgrid  models  are  supmior  to  the  mod- 
ds  widi  fixed  coefficients,  we  are  now  in  the  |nx>cess  of  evaluating  dynamic  subgrid  models  for 
compressible  flows.  Preliminary  results  show  good  agreemmit  with  the  results  of  Moin  et  al. 
(1S193).  More  results  of  this  study  will  be  be  n^iorted  in  die  near  future. 


3L  LES  OF  FLOWS  PAST  REARWARD  FACING  STEPS 


We  ue  now  getting  ready  to  samulate  more  complex  flows  such  as  flows  past  rearward  facing 
atq».  We  have  completed  preliminary  validation  studies  using  the  simple  eddy  viscosity  model 
(with  no  dynamic  evaluation)  and  have  draumstrated  the  ability  of  our  numerical  solver  to  repro¬ 
duce  result  ctmsistent  with  earlier  studies.  For  example,  Hg.5  shows  the  variation  of  reattache- 
meat  distance  (normalized  by  step  height)  as  a  function  of  Reynolds  number.  Also  shown  ate 
results  obtained  by  (Mher  researchers,  dearly,  our  LES  solver  is  in  good  agreement  with  earlier 
studies.  Figures  6a  and  6b  show,  respectively,  the  vortical  structures  downstream  of  the  step  for 
die  two  Reynolds  numbers.  As  Reynolds  number  iixrteases,  mote  complex  flow  patterns  are 
formed  as  expected. 

The  above  results  were  otMained  using  the  classical  eddy  viscosity  model.  These  earlio’  calcula¬ 
tions  were  carried  out  to  determine  the  accuracy  of  the  code  and  to  resolve  all  the  programming 
issues.  Therefore,  detailed  analysis  of  the  data  have  not  been  carried  out  Since  the  analysis  of 
subgtid  models  in  isotropic  turbulence  clearly  biggest  that  the  dynamic  models  are  superior,  we 
ate  now  in  the  process  of  including  the  dynamic  models  into  this  code.  It  is  expected  that  all 
future  calculations  in  this  configuration  will  be  carried  out  using  dynamic  models  (such  as  the 
dynamic  k-equation  and  dynamic  k-h  models)  and  for  relatively  high  Reynolds  numbers.  The 
exact  test  conditions  have  not  yet  been  finalized  siix:e  we  want  to  first  make  sure  that  there  is  suf¬ 
ficient  experimental  data  for  model  validation  in  such  flows. 

4.  PLANS  FOR  THE  NEXT  YEAR 

The  research  carried  out  in  the  first  year  focussed  on  simple  flows  such  as  isotropic  turbulence. 
The  methods  developed  for  analysis  of  subgrid  models  will  now  be  used  for  more  complex  flows. 
Thertfote,  from  now  on,  all  studies  will  focus  of  ctmiplex,  high  Reynolds  number  flows.  Two  test 
flows  have  been  chosen  for  subgtid  model  validation  studies.  The  first  configuration  is  the  rear¬ 
ward  fitcing  step  described  above.  The  second  configuration  is  a  co-axial  jet  shear  layer  with  and 
widiout  swirl..  This  type  of  highly  3D  swilling  flow  occurs  in  many  flow  situations  and  has  some 
interesting  features  associated  with  swirl  induced  mixing.  This  flow  is  also  sufficiently  complex 
and  there  is  stnne  experimental  data  for  comparison.  Both  these  configurations  will  be  studied 
using  the  same  numerical  solver  and  subgiid  moitels.  We  are  interested  in  determining  if  the  same 
type  of  subgtid  model  is  oqtaMe  of  simulating  accurately  these  two  different  types  of  flows. 
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[sosurfaccs  and  subgrid  hcliciiy  contours  obtained  during  LES  of  Taylor-Green 


Figure  2.  Contours  of  subgrid  helicity  computed 
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(b)  Compressible  scale-similarity  model 

Hgure  3.  Cross  correlation  analysis  of  compressible  subgiid  models  using  DNS  data  of  compress- 
fl>le,  (tecaying,  homogeneous  isotropic  turbulence  at  Re  =  10.  Subgrid  stress  correlation  obtained 
using  box  filter. 
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(b)  Compressible  scale-similarity  model 

Hgure  4.  Cross  correlation  analysis  of  compressible  subgrid  models  using  DNS  data  of  compress¬ 
ible,  decaying,  homogeneous  isotropic  turbulence  at  Re  s  10.  Subgrid  energy  transfer  correlation 
obtained  using  box  filter. 
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Figure  5.  Reattachmenl  length  computed  from  l^S  of  rearward  facing  step  flows  for  different 
Reynolds  number. 


(a)  Re  =  777 


(b)  Re  =  6992 

Hgure  6.  Instantaneous  spanwise  vorticity  contours  downstream  of  the  step  for  different  Reynolds  ^' 
munber  flows 
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EnENDED  ABSTRACT 

Ihe  coRBOt  hi  nhfrid  BMdds  B  otaeaciaDy  vateii  of  the  rinfiia  adify  vBoodQf 

■oddtotpropoaadbySmasoriBikyhilikB.  lUiitlBipteofihaMrioMltadiaiiantofaachMOddsthat 
hawbaaolwoaghtiBtotfiBaiioniDfaoeatMiidka.  toe»«yie.tact<iaiier  ofaaeiiy  ftoafteaBwaolvad 
to  Re  looived  acaks  hai  baca  ohaenrad  10  oocw  hi  dbact  MBarioal  riMritfioM  (MIS)  data  (PIOBDdli  at 
aL.  1991)  deaity  daBKnnndm  Rat  tha  ooaaaolvad  sshfrid  praoaaaa  onnt  ba  aoddad  by  a  poidy  AaB- 
paihe  aaachaidat  b  addhkm.  Re  "oonsmt”  (Ra  SoDaioriDBky  ooottaaO  appean  to  vify  Boa  Sow  to  flow 
ad,  Rea,  Re  aah|rid  model  leqahed  flsiher  Rie4DDing  for  a  tpectfc  flow.  Recem  ammpw  to  haprofve 
Ra  eddy  viaooaky  model  tavolvea  Re  adaptive  evatoatko  of  die  aahpid  coostaat  (Ganaane  at  aL.  1991) 
aad  ao  expBdt  amddhig  of  the  aobfrid  bachacager  process  (Chaauir,  1991).  b  apte  of  Reae  davelop- 
aMait,itlaBOtyetcleartfaiichahqdeiiiodelaaaeaofBckottorlaneeddyahDaiarinB(lES)ofhighRqr» 
BoldaBBmberfloiifs. 

The  tildecdve  of  ttda  paper  ia  to  evabate  the  eflMt  of  Re  fooB  of  Ra  choaea  aobfrid  aaodd  00  Re 
competed hwncaleaBeifytiaoaferprooeaabetweea the leaolvedaodaBaeaolvedacaieataLES.  Iheeveo- 
taal  foal  of  Ria  aaaaaRb  ia  to  davekfi  a  ad^rid  aaodel  Rat  wffl  periboa  adequate^  b  LES  of  bi|h  Rey- 
aoUa  maiRar  flowa.  1b  bveadfate  the  eflbct  of  Re  aobfrid  aaodd  00  Re  biBCKale  eoetiy  eaaafer  pro- 
ceaa,  me  b^  Re  atady  by  caaryief  oot  boR  ms  aod  LES  for  Ueeiical  fltm  oooRtkaa  aod  thao  aoriyz- 
htf  Re  eoeiiytaoaftrproceaaBaiafaoiathod  developed  aaceofly  by  PtwamdiM  at  aL  (1993).  Theyhave 
deoMoatmed  Rat  aah-frid  acak  ^OS)  eoerfy  tnoafer  qoaotitiea  cao  be  leadfly  exoacied  fiom  ms  data, 
b  taotaopmaiibk  flofw.  Re  coaffy  qwctram  £^(k)  of  Re  leaolved  acaka  at  waveomober  h  evolvea  by 

4eH») --2v**£^)+IXLH;,)+r(*ILeX  k<k,  (1) 

iihemvkReldBeaMticviacoaiy,i^kacet-offwaveoombcraqmBeiof  ReieaolvedaodaOb-fridacakt, 
TXkll^)  Npteaaott  eoeigy  ttaaafer  from  hoeiBCdom  whh  raatdved  acaka  ooly.  and  7*01  ii^)  rqeeaeatt 
btencRaMedRSGSmodea.  IfRefonoofiheqiectnlSGSed^vkooaityischoaeoas 
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rail.) 

iPFcw* 


k<K 


(2) 


fei  DNSl,  al  Mtt  ii  (1)  m  Ibkmm.  Itae  mt  hom  pnbkBs  with  this  approach  wha  LES  is 

tauUmi,  ii  LBS,  *•  kst  «■  hi  Bfatioa  (1)  is  aodeled  by  the  sriigrid  aodcL  1has.*ecMf|y 
IBMiv  IRMM  aar  te  cMitf  fwvoMaher  wfll  be  aflhciBd  by  the  fDOD  (or  choice)  of  the  Hhfihl  BOdd 
■Md  *•  alka  flf  iw  Mnolved  soiee. 

llii|aperwffliddtmihiipwiCMiarh»a^cfli»Biaf  theiesateten  WSspdUESaii^yi* 
e««dd^vlM0iiqrB0dBbaMh«*(Meshidiedc«llerbyChMBov(1991)aBdOenBnoeiaL(1991X  b 

■oM  to  to  ahbito  habtoB  hbetic  caeqy  ad  Mbfihl  hdidQf  (YtMhtowa.  1993)  wffl  be  evitaad. 
Mlher  esder  bbpld  aedeb  ae  aaldeied  eiietosi  to  touto^  high  KeyaoMi  assber  lows  toce. 
wIditoaitohleitIdieBetoiatoleatbacshi  toteMBaisaMsohedcacoatthieriyiHlreBfi  ewe 

«i#M— ■!«»  jy  fc  "y  V>  «li— ModdS  dMI 

ihtode  idbfrid  H***r**y  Otoe  toll  heUcity  ^weimd  Is  eacdy  wo  to  two-dbaakwl  swciBres) 

aqrbeaUeiBaooactotoeiyaaaiy  bitosBiifiidsadet.  These  esodels  wOl  be  aMbed  Mtaf  both 
tpecinl  SBd  physicsi  §faat  (flnise-rohae)  low  solven.  At  preseai.  to  netbod  of  adytii  is  sfanUar  to 
tot  inpoaed  by  Doandto  et  sL  (1993);  however.  «e  hope  to  dev^  a  eqiiiveleac  aneihod  of  adysis 
tot  ca  be  aed  dheedy  hi  to  phystol  qjwoe  ad  wfll  at  teqefae  hastomhn  the  dow  idd  too 
waeaBbcr  tpae>  This  is  pMlceliily  hoportat  ilooe  hi|h  Reyalds  lanbfr  dows  ot  pncdod  haiwt 
oopv  Ii  aapta  doada  (04^  nawitd  tobf  steps)  ad.  ba  M  will  at  be  possible  to  evdaie  to 
aetfy  iHMto  abi  Foaier  aomtaas. 

Ii  tos  piper,  we  coaider  dw  eoeipy  haste  properties  of  deo^fiof  hicaipettiible  ad  coefyresti- 
Me  isohcplc  hobileae.  b  the  doOowtap.  we  show  ad  coHyeie  pRUtoary  lealts  obtshifd  by  foor  df* 
teeaaedtods:  DKSabg  toooaiatdeashypieiid(HpectiddioiitoiofRofdto(1981).lb^abf  a 
caapeeasUe  dtoMielaa  code  (Meaa.  1992)  a  low  Mach  aator.  LES  aie|  •  coaiiwade  edify 
aodd  b  to  dedtovobae  aid  abs  i  iiniir  rtWrrrfTirf  a*  opdCx  The 

DNS  ito  LES  dsa  ae  Ofaataed  whh  64^  sod  32*  aid  ootoL  tesoecdvdy.  fo  dftecnoes  b 

aaafcd  asobda  ad  aetods.  to  paasetea  see  choea  a  shadsa  stadadcalfy  to  saae  dow  b  an 
caes.  MoadetdboftoamerlalaethodsadtosobttMaoddswfflbeiivabtodBdpspcr. 


Bcthabf  toa  a  botapb  Oaessia  andoa  ddd  whh  a  padded  tahbl  euagy  spectnmi,  to 
bjibnifyaimlr  dcM  b  allowed  a  evehe  aadl  a  *aaliadc”  eelf-stadbr  stab  b  developed  (a  b  Yeoof  and 
Basaaa  1991).  Resohi  toa  a  6^  cakabdoii  tow  a  powerbw  decay  of  eaeny,  with  eqaaea  131 
cceebtea  whh  |rid  peaeraad  tatobace  dea.  ad  saaU-ecab  adveraUfy  iDmtrsted  by  to  colbpse  of 
bglHahvaaBtor  epecha  a  dUteca  tbas  aader  tobDOforov  scalbi.  a  town  b  HfURs  lad.  bb 
baeatbi  a  nob  dai  both  to  DNS  seshhs  agiee  loaaikafaly  wdl  and.  fbihennoie.  dat  to  LES  dab  abo 
tosa  vcqr  pood  atacaaa  wib  to  DNS  dab. 


The  dow  teUs  a  vabes  daws  wea  welyad  Here,  to  lesrii  obtained  a  a  dae  i«9X>  (vdiich 
caaBtoBd*aaoaad203bldaadd)ManovatlBa)waadeciedtotatoranafysb.  Bydibdae.to 
tatobnoa  b  wel  devdoped.  aad  to  Tbgdorecik  ReraUs  aaaba  (baaed  on  to  r jnx  vdochy  and  Thy- 
tarafcroicddbaTwabaafylKiiowfydecaabiedbdaa).  FipHa2a>2dabowstotopaorto 
■oaaHaed  aaeqy  ad  dbefyatba  tpeam  a  to  choacn  tba  fbr  to  ten  dUfaea  caaa.  The  DNS  dab 
bw  fctiifc  Bt— SMtoM  MB  —a  — rfiint  wmA  mtut.  ftlfc  —A  fwrtFMy  mmO  The  pidCS  Of 

awfy  and  dbhfyadon  apecaa  ae  aqnnied  by  a  toior  of  two.  The  LES  daa  abo  afaa  aaeonhUy  wib 
to  ONS  laaahs.  howBvw.  to  peak  vahie  of  both  to  eneror  and  dbsipaioo  b  sdihdy  higher. 

flpae  3, 4  ad  S  tow  to  eddy  dsooshy  (Eqoacha  (2)X  sabgrid  nanste,  7*(tli^)  end  icaolved 
ttanateriifcli^Xnapecdvdy.  to  dtffcm  cotoff  wavennBd)ers.4«.  AhhonfhtoiaiteofKaksavailalde 
b  oa  64*  cikdaiow  b  tahed.  newtoka,  sane  hnereetbg  oonctasioiis  aay  be  dawn,  to  be 
oendnnedby  ISS’sesahitocwillbeobtataedby  tothneoftoSyntoostaa.  Rendtttotoloweaait- 
off  chosen  (1^  ai^  aa  aoaewha  eoade.  becana  thqr  ooRtqxad  to  LES  with  a  vciy  ooatse  dher.  and 
becaaa  of  etadstkd  varihbilify  aesnriaiwl  wib  to  lesoivedtob  eaaiy  iftanm  (E*’(k))  a  low 
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wMCBoaBbcn.  For  hi||wr  cniolb,  ipeciral  eddy  viscosity  ipcneies  teply  as  the  aaoB  wavemanber  is 
ifpraadMd  (Ftinre  3»c).  Note  that  ffar  die  LES  case  (Rfote  3c).  tecaioffwaweBnBiberhe>lS  is  actndly 
dtt  aHKiMai  ■»•«—«"«*«**  hi  tte  lesolved  field.  At  «BaD  k/k,,  tte  spectral  eddy  vlsoo^  Is  amative, 
dwwtagaahdtanatllBlialloaoftheeddly  visooshyooaoqithaelf.  These  fBatares  are  ocosliieat  with  Oe 
SMahi  of  Doawaddd  et  aL  (1993)^  who  studied  a  Thylor-Onea  vortex  nd  peeseaied  the  variadoa  of 
Va9l3;t)>hHCaoc*e«Mftn  (7XIII^)aod  r*(kli^))asafiaKtiaBorh/i^. 

Rfun  4  Aows  tel  the  eOeci  of  sub-frid  scales  is  aalBly  to  enact  caeriy  fine  drt  hK|e  scales,  but 
then  Is  also  a  aoDHt^figMe  badkscatter  (r*(b  li^  >  0)  eOtet  acthtf  OB  fee  kcfett  scales.  We  pIsB  to  rtady 
fitherthliclhirthyeayacldyaeoddhnthehafhiramrprocessasdDBeeartertyCltasaov(1991).  Flgare 
3  shows  that  the  aatnre  of  laieiactloasaBMat  the  rcecdved  scales  thCBrtcives  It  priaiarilyfarwrdcaicad- 
lag:  wMh  eae^DT  lost  (7UII«)<Q>  St  low  h/Aif  sad  caei|ytahi(rOlli^)>0)  as  h/l^  approaches  aahy. 

These  lesalts  wifi  be  Ihrtbervslhisied  far  the  fiaalpeper  by  canyhn  out  128*  DNS  shanlarinni  As 
aotsd  above,  aubgrid  aeodels,  each  as  die  eddy  visooshy  aeodd  whh  adapdve  evahniiaa  of  the  ooastam 
COeesusw  et  aL.  1991)  sad  the  oao-e^iailoa  aad  two^^ndoa  aaodela,  will  be  isnpkmenitid  hi  bodi  the 
fialie>vofaBBe  oodes  tsid  the  effect  of  te  aabgiid  modd  flo  dw  cBcrgy  asBsftr  prooeu  win  be  tavesdgaied 
la  detail  adag  a  aeries  of  LES  at  vaifcMS  grid  lesohakxi.  The  effect  of  locrcasiBS  the  Rcyaolds  otaaber 
wifiaiaobelavestlgaied.  The  results  of  these  studies  should  abed  Ugbtoa  the  defideades  sad  strcagtbs  of 
fee  sabgrid  awdeh  aad  provide  a  dfaeodoB  far  hnproviag  the  saodelt  so  that  high  ReyuoUt  auadier  fiows 
csB  be  aiiaalsied. 
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l.lNIRODUCTiON 

Tine  (fynamic  idbgiid  teak  (SGS)  aaodding  technkyic.  iatroduced  by  Getmano  et  al.  (1S>91). 
has  been  snooessOiUy  NVUed  to  various  types  at  flow  Odds  (see  Moin  et  al..  1994  for  a  leoent 
levlew).  There  mt  three  iopoitant  features  of  this  modeling  approach.  Rrst,  the  modd  constant. 
the  so-called  Smagoriasl^'s  constare,  is  not  prescribed  a  priori,  father,  it  is  determined  as  a  part 
of  the  solution.  lUs  sppeoach  removed  one  of  the  mi^  probtems  assodated  wife  fee 
SBUga>li>sky*s  eddy  vlaoosity  modd  which  was  the  determl  nation  (and  One-tunios)  of  the 
constant  for  dUfetent  flows.  Seemd,  the  dynamic  modd  diows  fee  correct  asynqjtotic  behavior 
near  wans  and  third,  as  a  resuh  of  the  dynamic  evahurioo.  the  oonsuau  can  become  negative  in 
certain  felons  of  the  flow  Odd  and  feus.  Nspeats  to  have  the  c^MbiUty  to  mimic  backscatter  of 
energy  ftom  the  solved  scales  fo  the  resolved  scales.  Ihis  last  feature  is  particolariy  attractive 
dnee  direct  numerical  simulation  (DNS)  data  has  shown  that  backscaiter  effect  can  domioate 
over  a  signiflcaat  portioo  of  the  grid  potans  In  the  flow  ilidd  (e.g..  Tloindli  et  al.,  1991). 
However,  large-eddy  sjmulatkmsQLES)  using  this  dynamic  inodd  has  dtown  feat  fee  bacitscaner 
effect  of  OetmanoY  modd  era  be  excessive  raid  can  cause  numerical  instability.  Some  recent 
efforts  (Moin  et  al.,  1994;  Ffomdli  and  Liu.  1994)  were  directed  particularly  to  develop  better 
methods  to  evaioate  the  cofutonr.  However,  these  studies  atfll  assume  that  dte  form  of  fee  eddy 
viaoosity.  as  proposed  by  Smagotlnsl^.  is  sdn  appropriate  and  valid  for  fee  entire  range  of  flows. 
Reoeai  IMS  studies  of  decaying  isotropic  natulenoe  suggest  feat  fee  assumption  of  fee  local 
equUibriam  between  fee  SGS  energy  prodnetiao  and  dissipation  rate  (whidi  was  used  to  derive 
fee  Smagorinsky  modd)  may  not  be  satisfied  over  a  large  portioo  of  fee  grid  poiats.  Ibis  lack  of 
local  equilibrium  can  be  significant  wba  the  flow  Adds  at  different  grid  resohitioos  are  rdated 
10  eadi  other,  as  is  the  case  in  fee  dynamic  modd. 

la  tUs  paper,  these  Issues  are  addressed  by  oonsMetliig  two  new  subgrid  modds  that  do  not 
require  locd  equilibttom  between  the  SGS  energy  production  and  dissipation  rate.  Iheae  modds 
have  been  devdoped  using  the  dynamic  moddng  approach  and  iherefere.  there  are  no  coRSfunu 
that  have  to  be  prescribed  a  priori,  Ibe  first  modd  feat  is  studied  is  a  dynamic  versioo  of  fee 
KdiBMigotQfvY  scaling  exptearion  fbr  eddy  viscosity  nd  the  second  modd  is  a  ooe-equdioo 
reodd  fer  SOS  kinetic  energy. 

IMag  LES  of  decayirtg  isoiropk  mrbuleaoe  at  various  Remolds  annfeer.  the  behavior  of 
these  new  tfyaniic  aiodds  have  been  compared  to  fee  pretUcti^  of  the  dynamic  modd  based . 
on  fee  Smagotinsky's  eddy  viaoodty  oooc^  An  evetmial  goal  of  this  research  is  to  determine 


tfie  ftm  of  a  dynamic  lubgrid  modd  that  wlU  provide  realistic  results  over  a  wide  range  of 
EeynoMs  number  and  grid  resohitioo. 

Tire  numerical  simulatioos  were  carried  om  using  a  llniteKliffiereoce  code  dut  is  second-order 
accuracy  in  time  and  fUlb-order  (die  convective  terms)  and  sixth-order  (the  visoous  terms) 
aocmacy  in  ^reoe  using  upwind-biased  diffemoes  (Rai  nd  Moin,  1991).  Time-accurate 
iohidoiis  of  the  Incompressible  Navier-Stokes  equations  are  obtained  by  the  artificial 
compressibility  approach  (Rogers  <r  u/..  1991)  which  requires  subiteration  in  pseudotime  to  get 
dre  diveqienoe-flm  flow  flekL  Earlter  (Menoo  and  Yeung.  1994).  this  code  was  validated  by 
onying  out  IWS  of  decaying  isotropic  furbulenoe  and  cooqnring  the  resulting  statistics  with  the 
predictioos  of  a  wdl  known  psuedo  spectral  code  (Rogallo.  1981). 


Z  DYNAMIC  SUBGRID  MODELING 

In  the  Hollowing,  the  ftnns  of  die  dynamic  subgrid  models  studied  here  are  briefly  described. 
More  details  of  the  moddlng  approach  will  be  given  in  the  flnal  paper. 

2.1  Dynamic  Smagorinsky  Model 

The  primary  requhenrent  in  LES  is  the  modeUag  of  the  SGS  stresses  that  result  frran  the 
spatial  flltering  of  the  instantaneous  Inoompressible  Navier-Stdto  equations.  If  ao_^appropriate 
filler  is  eaqiloyed  (here,  a  top  hat  fliter  with  a  filter  width  A  is  used,  where.  Als  the  grid 
spacing),  dren  the  true  SGS  stress  tensor  is:  This  stress  tensor  must  be  modeled 

to  close  die  LES  equations  of  motioiL  The  Smagorinsky  eddy  viscosity  dosure  for  the  SGS  stress 
tensor  is  of  the  form: 

•  Vt«CA*P|  (1) 

ndrereC  is  the  model  coefficient,  5^  is  the  resolved  scale  strain  rate  tensor 


and  jJ]  B  (2^^)’^.  Here,  the  ovnbar  on  the  flow  vreiables  indicates  die  result  of  flltering  at 

the  grid  scale  A.  In  die  dynamic  moddlng  approach,  a  mathematical  identity  between  the 
serenes  r^ved  m  the  grid  scale  Alter  A  and  a  test  filter  A  (typically  diosen  to  be  twice  ctf  the 
grid  fliter  A)  is  used  lodeienninethe  modd  coeffidemCtt  a  port  of  the  simulation.  Tims,  if  the 
application  of  die  test  filter  on  any  wUble  4  is  denoted  by  ^  or  then  it  can  be  shown  that 

I*  ■ 


(3) 


Here,  («  Is  deflned  uslBg  the  test  filter.  Using  file  assunqxioD  of  fiie  self- 

stanUaiity  of  the  sifitgiid  stress,  7^  can  be  modded  in  the  same  way  as  resulting  in  die 
following  exiiression: 


“  2  v,  ^ 


Vj  *  Cd? 


(4) 


SidMtituting  (1)  and  (4)  into  (3),  one  can  obtain  an  equation  for  C 


(5) 


ndiere 


(6) 


Equadoo  (5)  is  a  set  of  indeoendent  equations  for  one  unknown  C.  To  minimize  the  error 
ftom  solving  ttiis  ovetdemrmined  system.  Lilly  (1992)  proposed  a  least  square  method  «4iich 
yidds 


c^IMl 


(7) 


Past  studies  (e.g.,  Zang  et  al.,  1992;  Yang  and  Peiziger,  1993)  using  this  qiproadi  have  shown 
fiiat  die  vdue  of  C  obtained  ftom  Equation  (7)  can  vary  widdy  in  die  flow  £Md.  Tliis  can  cause 
numerical  instability.  To  idieve  diis  pcdUem,  qiatial  averaging  is  typically  performed  for  both 
die  numerator  and  die  denominator  on  the  RHS  (7)  (Piomdli.  1993).  Usually,  diis  averaging  is 
done  ody  in  the  difectkms  of  flow  homogeneity  (e.g.,  Moin  et  al,  1991).  In  die  present  study  of 
homogeneous  isotropic  turbulence,  averaging  is  implemented  over  the  entire  oomputationd 
domain,  hence  C  is  a  flinction  of  time  only. 


2.2  Dynamic  Kolmogorov  Scaling  Model 

RecenUy,  Wong  and  Lilly  (1994)  reevaluated  the  KolmogorovOs  scaling  expresston  for  die 
eddy  visoodor.  Hie  representative  wiables  characterizing  SOS  fluctuations  are  the  SOS  energy 
disaipalion  cate  e  and  die  SOS  lengdi  scale  rriiich  can  be  approximated  by  the  grid  interval  A. 
Usiiv  dtanensknai  arguments,  the  subgrid  eddy  viscosity  can  be  written  as, 

Vt  *  C,A^  (8) 

Here,  the modd  ooeffident  C,  has  the dlmensioo  rdated  to  e  Qaecisdy,  "E^).  Ihis  rinqde 
ciqiession  ibr  can  be  robust  because  it  does  not  enqfoy  the  assumption  ttf  local  equilibrium 
between  SOS  energy  production  rate  and  dissipatioo  tare  nhldi  was  adopted  to  derive  die 
Smagoiinslgr  modd  (and  udiich  is  iiqilidt  in  the  (tynamic  modd  described  in  Section  2.1). 
Another  advantage  of  tUs  modd  is  that,  due  to  die  sinqfidty  of  die  model,  die  total 


OQoqMtttknal  time  can  be  algnlflcanfly  mhiced.  However,  this  modd  can  be  valid  only  by 
enqiloylng  a  dynamic  modding  method  to  determine  the  modd  ceffident  Otherwise,  there  may 
be  no  way  to  prescribe  ttiis  non-dlnmensiooless  ooefllclent  by  a  fixed  value.  The  resulting 
equatioa  Ibr  C.  is  eaqpressed  as  follows. 


2.3  Dynamic  t-equatkm  Kfodd 


A  higher  order  modd  is  also  oonsidred  in  this  study.  Ihis  is  a  one  equadcm  modd  for  die 
SOS  Unetic  energy  k  (m  (uf-u^^}/2).  This  modd  was  recently  evaluated  by  Menem  and 

Yeuog  (19SM)  through  n  priori  tests  using  IMS  of  decaying  and  forced  isotropic  turbuloice.  In 
this  mo^  an  evotution  equation  for  H  can  be  written  in  fire  following  form 


where 

--2VjS^-h^8^k  ,  Vj-Qk'^^A  (11) 

3 


Also,  one  can  model  the  (Ussipatkm  rate  e 


(atV  (dui/d!Xj)  ) using ir,  as 


(12) 


Menon  and  Yeung  (1994)  chose  the  modd  ooeffidents,  based  on  an  earlier  study  (Yoshizawa, 
1993),  to  be  C„«0.0854  and  C(«>0.916.  However,  these  values  were  found  to  cause  excessive 
diss^wiitm  of  large-scale  energy  for  isotropic  turbulence.  In  die  presem  study,  we  apply  the 
dynamic  modding  tedmlque  to  obtain  qiproptiate  values  of  the  ooeffidents.  To  inqtlement  die 
dynamic  approach,  the  sobgdd  Unedc  energy  JT  at  the  test  filter  levd  is  required.  This  is  obtained 
by  using  die  trace  of  (3),  K  s  I^/ 24  ft.  Thus,  using  the  procedure  described  earlier,  equadems 
for  bodi  Q  and  C,  can  be  derived: 

-A(ft'^  J,))  (13) 


(14) 


Here,  equation  (14)  Is  a  acdar  equation  fior  a  scalar  unknown  and  tluis,  we  can  obtain  tiie  exact 
vdue  Ibr  C,  wWioitt  api^ylag  tte  least  square  medxxL 


OnemaytliInkthattlieinetttofttiBLESapproecb  lies  in  the  slnq^city  of  the  model  used  for 
tlieaiiiiulations  and  that,  this  merit  may  be  lost  by  using  the  higher-order  modds  such  as  the  one 
■equation  modd.  However,  die  increased  cost  of  this  type  of  modd  is  compensated  by  the 
improvement  in  the  accuracy  resulting  by  not  assuming  local  balance  between  the  SGS  energy 
production  and  dias^ation  me.  Aoooniing  to  our  n  priori  test  using  rdativdy  higb  resolutiofl 
DNS  (usinf  128x128x128  and  64x64x64  grid  resohitfons),  over  a  large  portion  of  die  grid  points 
in  iaoiropic  turbulent  flow,  local  equillbrinm  assumption  is  vioiaied  Ihis  implies  that  modds 
flmt  do  not  lequiie  local  equlibrium  such  as  the  one-equation  modd  has  ttie  poientid  to  produce 
better  results  tiian  die  Smagorinsky-^pe  modd.  Ibe  earlier  studies  using  die  one-equation  model 
widi  fixed  values  of  the  coeffidents  (Menon  and  Yeung,  1994)  Improved  the  results  vriien 
compared  to  the  Smagorinsky’s  modd  with  fixed  value  of  coefflcieat.  However,  the  results  were 
poorer  dm  the  results  obtained  using  the  tfynamic  Smagorinsky*s  modd.  Analysis  of  the 
stanidatioo  data  dwwed  dut  this  was  (hie  to  a  poor  predictiott  of  bodi  the  production  and 
tflssipation  tenns  in  (10)  using  fixed  ooefBdems.  Eqie^y,  the  prediction  cf  the  dissipatioo 
terms  was  very  poor.  In  the  presem  study,  die  use  of  the  dynamic  procedure  to  evaluate  die 
ooeffldents  results  in  a  mudi  better  predictioo  (rfboth  the  production  and  dissipatioo  terms  in  the 
k-equation  modd.  Hds  in  turn  inqxoves  significantly  die  results  of  the  LES  udng  the  dynamic  k- 
equation  subgrid  modd. 

3.  PRELIMINARY  RESULTS  AND  DISCUSSION 

To  evaluate  the  three  tUflerent  dynamic  SGS  modds  (described  in  Sec.  2),  LES  of  decaying 
hoanogeneous  isotropic  turbulence  were  conducted.  Stuting  from  an  initially  divergent  free 
vdodty  Add  with  a  prescribed  energy  qiectrum,  the  flow  is  allowed  to  devt^  into  realistic 
decayi^  turbulence.  Ihe  flow  under  consideratioo  is  modded  in  a  cubic  box  with  periodic 
boundary  conditions,  and  two  grid  resolutions,  32x32x32  and  16x16x16,  were  employed  for  the 
simulations.  Die  stanulations  were  performed  for  Oree  difrerem  initiai  Taylor  Reynolds  mindiers 
Re;L,*30,  100  and  1000.  In  ttie  lowing,  prdiminary  results  are  discussed  to  highlight  the 

behavior  of  die  new  dynamic  subgrid  models.  Die  results  discussed  here  are  obtained  by  using 
the  flow  Add  at  an  instant  in  time  vriiere  die  flow  Add  has  rdaxed  to  a  realisiic  decaying 
tmbulence  and  die  Reynolds  number  of  the  flow  Add  is  decieasiitg  very  slowly. 

To  evaluate  the  sdf  consistency  of  the  dynamic  models,  the  LES  results  obtained  on  the 
32x32x32  grid  resolution  were  compared  to  the  LES  results  obtained  on  the  16x16x16  grid.  To 
ensure  that  both  grid  resolutioo  simulatioos  are  being  performed  using  nearly  ktemical  initial 
conditioos,  the  initial  flow  Add  for  die  16x16x16  grid  atanulation  is  obtained  by  Altering  die 
32x32x32  grid  resolved  flow  field.  Die  resulting  flow  Adds  from  these  two  diffierem  resolution 
simulatioos  can  be  idaied  usta^  die  madiematicd  identity,  (3).  By  dds  approach,  die  modeled 
quantity  can  be  obtained  from  the  two  data  aets  which  dxwld  be  ideotical  or  at  least 

bfgbly  coirdated  if  diesiAgrid  inodd  performed  correctly  at  die  two  grid  levds.  To  tyiantiiy  die 
behavior  of  die  subgrid  model,  cocrdation  coefficieitfs  (defined  in  die  usual  manner)  are 

oonqmted  usiitg  the  anisotropic  parts  of  esqxessioiis  for  Die  variation  cf  die  averaged 

condationooefBcieiasofthese  three  anisotropic  componentt  are  shown  in  Fig.  1  as  a  ftinction  of 
the  Taylor  R^ndds  number  at  the  instant  of  the  oom^soo.  Die  coRddloo  coefficient  is  very 
high  for  an  the  models  over  the  entire  range  cd  Taylor  Reynolds  numbers.  However,  the 


ooRdalioa  coefficient  for  foe  dymmic  Snugocinsky  modd  decreases  with  increase  in  the 
Reynolds  numbers.  Altboo^  the  ootrdatioo  coefficient  for  the  dynamic  Kdmogorov  scaling 
modd  is  always  dian  the  others,  it  was  found  using  qwctral  qiace  analysis  diat  this  model 
has  some  problems  in  pretUctiiig  the  dissiiMlioo  rate  of  the  SGS  turtmlent  energy.  This  modd  is 
more  diss^Mtive  in  the  earlier  transitional  stage  and  less  dissave  In  the  folly  devdoped 
ttibulent  stage.  This  problem  m^  be  caused  by  the  lack  of  direct  modding  d  the  dissipation 
nm  svhen  die  KolmogorovOs  scaling  for  the  edify  viscosify  is  used.  In  any  evem,  the  results 
shown  In  Fig.  1  dearly  suggest  diat  the  modds  that  do  not  make  the  assmqition  d  local 
oqnlllbilttm  between  the  production  and  diss^on  rate  produce  better  resutts  than  the  dynamic 
Sma^xlnsiys  modd.  This  is  Impcntant  since  the  eventual  goal  of  LES  medioddogy  is  to 
devdop  su^tid  modds  that  wOl  aUow  ttanuladon  of  high  Reynolds  number  flows  using 
fdadvely  coarse  grids  (grid  icsoludoo  lestrictions  are  typically  imposed  due  to  computer 
leaouroe  limitations). 

Fig.  2a  dnws  the  variadon  in  die  dynamically  evaluated  constants  with  dme  during  die 
ahmdadon  for  Rcj^^lOO  and  Fig.  7b  diows  the  variation  of  the  modd  constant  widi  die  Taylor 

Reynolds  numbers.  Obviously,  die  modd  oodBdems  go  through  changes  in  die  earlier 
transidooai  stage.  However,  after  some  time,  dl  ooefOdenis  reach  an  asymptotic  state  with  the 
excqition  of  C,  in  Kolmogorov’s  scaling  expression  which  keqis  decreasing  as  the  kinetic 
energy  decreases.  Also,  the  vdues  of  the  coeflldents  at  this  asynqnotic  state  are  dmost 
tinVpMMVau  on  Reynolds  number  except  for  C,  in  k-equadon  modd  because  E .  and  hence  Q 
vdddi  is  generated  ftom  direct  modding  of  e,  is  very  sensitive  to  the  grid  lesohition  and 
R^ndds  number.  The  resulting  ape  dmflf  to  the  vdues  obtained  in  earlier  studies. 

For  exaaqile,  in  the  present  study,  die  tfynamic  Smagoiinsky  modd  predicts  that  the 
Smagorinsky  modd  ooefSdent  C,  (which  is  the  square-root  of  dynamic  modd  coefficient  O 
Should  be  about  0.16S.  lUs  is  very  dose  to  the  vdue  of  0.17  suggested  by  Lilly  (1966)  for 
homogeneous  isotropic  turbulenoe  with  cutoff  in  the  inertid  subrange.  It  is  noteworthy  that 
detetming  the  modd  ooefficient  using  die  dynamic  Smagorinsky  modd  may  have  some 
limitations.  This  modd  predicts  its  modd  coefficient  to  be  nc^ve  for  a  long  period  of  time  in 
the  earlia  transitiooal  stage.  Even  this  period  of  time  Increases  with  increasing  Initial  Reynolds 
number.  Obviously  it  leads  to  numerical  InstdtiUty.  To  prevem  diis  problem  from  happening,  we 
inqiose  die  constraints  diat  C  should  be  dways  lar^  dian  0.01.  This  kind  of  numerical 
Instability  was  not  brought  about  by  die  other  modds. 

Fig  3a  shows  contours  of  die  SGS  kinetic  energy  on  32x32x32  grid  dhecdy  obtained  by  LES 
using  die  dynamic  k-equatkm  modd.  This  result  Is  cotnpired  to  die  predicdon  by  the  LES  using 
the  same  modd  but  diffornt  resolution.  16x16x16  grid,  (Fig.  6b)  at  an  artitrary  (but  same)  slice 
of  the  3D  field.  The  oonqiaiison  diows  that  there  is  significaiit  similaiity  between  two  LES 
results  of  dUferem  grid  widths  and  dnis  the  sdf  consistency  of  die  dyruunic  k-equation 

modd.  The  peak  vdues  of  the  coarser  grid  (16x16x16)  results  are  approximatdy  twice  of  those 
ftom  die  finer  grid  (32x32x32).  This  is  reaaondile  because  the  coarser  grid  nhidi  has  die  lower 
cutoff  wave  number  diould  contdn  die  larger  SGS  kinetic  energy  indde  of  its  subgrid  regions  for 
the  same  Reynolds  number. 


4.  FUTURE  WORK  TO  BE  INCLDKD  IN  THE  FIN  AL  PAPER 


Ibe  ffesuks  obt^ned  so  te  ri»w  that  the  dynamic  modds  ttiat  do  not  assume  local 
eqiiDiMum  between  the  SOS  energy  productioo  and  dissiiMiion  r«e  perform  significandy  better 
than  the  dynamic  modd  based  on  the  dassical  Smagoiinsky's  eddy  viscosity  model.  Funher 
aimMlaiions  are  underway  to  oonfinn  this  result  using  finv  grid  resolution  and  for  higher 
Rqmolds  nunrijers.  Gompailsoo  of  die  LES  results  with  the  NSIS  results  for  iqxesentmive 
Cmattvdy  low  Reynokh  numbers)  cases  will  also  be  carried  out  FfauUy,  these  models  will  be 
cwaiiialed  for  dlflbrent  flow  cases.  For  example,  die  Taykr-Green  vortex  flow,  which  has  been 
■aed  by  Domaradxki  et  al.  (1993),  is  ooosidesed  a  g^  test  problem.  Ibis  Is,  due  to  flow 
qnninetry,  a  flow  at  a  rdadv^  high  Reynolds  number  can  be  simulated  using  rdadvely  coarse 
feaoiution.  The  results  of  these  studies  will  be  described  in  more  details  in  die  final  paper. 
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